Majorization in Quantum Adiabatic Algorithms 
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We study the Majorization arrow in a big class of quantum adiabatic algorithms. In a quantum 
adiabatic algorithm, the ground state of the Hamiltonian is a guide state around which the actual 
state evolves. We prove that for any algorithm of this class, step-by-step majorization of the guide 
state holds perfectly. We also show that step-by-step majorization of the actual state appears if the 
running time becomes longer and longer. This supports the empirical viewpoint that step-by-step 
majorization seems to appear universally in quantum adiabatic algorithms. On the other hand, the 
performance of these algorithms discussed in this paper can all be estimated, which is exponential 
in the problem size. This can be looked as a strong evidence that step-by-step majorization is not 
a sufficient condition for efficiency. 
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I. INTRODUCTION 

In the past two decades, quantum computation has at- 
tracted a great deal of attention, because it was demon- 
strated that the performance of quantum computers ex- 
ceeds that of classical computers for some computational 
tasks. Among quantum algorithms proposed so far Shor's 
factorization algorithm [jj and Grover's search algorithm 
are two famous examples. However, it is believed that 
the design of quantum algorithms seems to be very diffi- 
cult . Thus studying the characters of quantum algo- 
rithms is a very important question. It has been observed 
that majorization theory seems to play an important role 
in the efficiency of quantum algorithms [j, [a, [6j • For ex- 
ample, in many quantum algorithms the initial state of 
the system is an equal superposition state and the final 
state before measurement is some computational basis 
state. In the process of computation, the probability dis- 
tribution associated to the state of the system in the 
computational basis is step-by-step majorized until it is 
maximally ordered. In 6] by carrying out a systematic 
analysis of a wide variety of quantum algorithms from the 
majorization theory point of view, the authors concluded 
that step-by-step majorization is found in the known in- 
stances of fast and efficient algorithms, such as quan- 
tum fourier transform, Grover's algorithm, the hidden 
affine function problem. By contrary, in j(| the authors 
also offered a quantum algorithm that doesn't show step- 
by-step majorization, which doesn't provide any compu- 
tational speed-up. These facts show that step-by-step 
majorization seems to be necessary for the efficiency of 
quantum algorithms. 

In Q and [U, the authors also studied majorization in 
quantum adiabatic algorithms. Quantum adiabatic com- 
putation is a novel paradigm for the design of quantum 
algorithms. In a quantum adiabatic algorithm, the evolu- 



tion of the quantum register is governed by a hamiltonian 
that varies continuously and slowly. If the initial state 
of the system is the ground state of the initial hamil- 
tonian, the state of the system at any moment in the 
whole process of computation will differ from the ground 
state of the hamiltonian at that moment by a negligible 
amount. Thus, in a quantum adiabatic algorithm the 
ground state of the hamiltonian is a guide, and the ac- 
tual state of the system is always around this guide. (In 
this paper we call ground states in quantum adiabatic 
algorithms guide states.) In Q and [jfl, through numer- 
ical simulations to several special cases the authors got 
an empirical conclusion that quantum algorithms based 
on adiabatic evolution naturally show step-by-step ma- 
jorization provided that the Hamiltonians and the ini- 
tial state are chosen with sufficient symmetry and the 
evolution is slow enough. In this paper, we will study 
the majorization arrow in a big class of quantum adia- 
batic algorithms. We prove that in any algorithm of this 
class step-by-step majorization of the guide state holds 
perfectly. For actual state, we show that step-by-step 
majorization appears as long as the running time of this 
algorithm becomes longer and longer. Based on these 
results, we analyze the relation between the efficiency of 
quantum algorithms and majorization. 

The rest of the paper is organized as follows. In Sec. 
II we briefly review quantum adiabatic computation and 
majorization theory. In Sec. Ill we prove that step-by- 
step majorization of the guide state holds. In Sec. IV 
we discuss step-by-step majorization of the actual state. 
Finally, in Sec. V we discuss the relation between the 
efficiency of quantum algorithms and majorization, and 
summarize our conclusions. 



II. QUANTUM ADIABATIC ALGORITHMS 
AND MAJORIZATION 
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For convenience of the readers, let us recall quantum 
adiabatic computation and majorization theory. 

Quantum adiabatic computation, proposed by Farhi 
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0, is based on quantum adiabatic evolution. Suppose 
the state of a quantum system is |^>(t))(0 < t < T), 
which evolves according to the Schrodinger equation 



H(s) = {l-s)Ho + sH u 



(9) 



i-\m)=H(t)\m), 



(i) 



where 



where H(t) is the Hamiltonian of the system. Suppose 
Hq = H(0) and Hi — H(T) are the initial and the final 
Hamiltonians of the system. Then we let the hamiltonian 
of the system vary from Hq to H\ slowly along some path. 
For example, an interpolation path is one choice, 



H(t) = f(t)H + g(t)H 1 , 



(2) 



where f(t) and g{t) are continuous functions with /(0) = 
g(T) = 1 and f(T) = g(0) = (T is the running time of 
the evolution). Let \Eq, t) and \E\, t) be the ground state 
and the first excited state of the Hamiltonian at time t, 
and let E (t) and Ei(t) be the corresponding eigenvalues. 
The quantum adiabatic theorem Q shows that we have 



|(£ ,TW))| 2 >1 



provided that 



D„ 



g 2 



<e, < e< 1, 



(3) 



(4) 



where g m in is the minimum gap between Eq (t) and E\ (t) 



g„ 



= min [Ex(t) - E (t)}, 

0<t<T 



(5) 



and D max is a measurement of the evolving rate of the 
Hamiltonian 

dH 

D max = max \(E 1: t\ — \E ,t)\. (6) 

o<t<T at 

Quantum adiabatic computation is a novel paradigm 
for the design of quantum algorithms. For example, 
Quantum search algorithm proposed by Grover [2] has 
been implemented by quantum adiabatic computation in 
0- Recently, the new paradigm for quantum computa- 
tion has been used to try to solve some other interesting 
and imp ortant problems, such as Deutsch-Jozsa problem 
[Tol llll Il3 , hidden subgroup problem Hll , 3SAT prob- 
lem 0,13) traveling salesman problem [15| and Hilbert's 
tenth problem [l6| . 

Let's look a big class of typical quantum adiabatic al- 
gorithms, which is what we will discuss in this work. Sup- 
pose / : {0, 1}™ — > R is a function that is bounded by a 
polynomial of n. Let Hq and Hi be the initial and the 
final hamiltonians of a quantum adiabatic evolution with 
a linear path H(t). Concretely, 



Hq = I 



N 



\a)(a\, 



(7) 



(8) 



1 N 

* i— 1 



(10) 



and s = s(t) a continuous increasing function with s(0) = 
and s(T) = 1 (T is the running time of the quantum 
adiabatic evolution). According to quantum adiabatic 
theorem, this class of algorithms can be used to minimize 
the function f(i),i — 1, 2, N (we suppose there is only 
one minimum). The quantum adiabatic algorithms for 
search problem in 01, hidden subgroup problem in |13| . 
3SAT problem in |l4| and traveling salesman problem in 
[THj l belong to this class. 

Now let's turn to the majorization theory. Majoriza- 
tion is an ordering on N-dimensional real vectors. Sup- 
pose x = (xi,x 2 , ...,x N ) and y = (yi,y 2 , are two 
N-dimensional vectors. If x is majorized by y, y is more 
disordered than another. To be concrete, let mean 
x re-ordered so the components are in decreasing order. 

We say x is majorized by y, namely x -< y, provided 
k k ' N N 

E 4 < E v\ for * = 1.2, ...,N- 1 and £ x\ = £ y\. 

i—l i—1 i—1 i—1 

It has been proven that majorization is at the heart of 
the solution of a large number of quantum information 
problems. For example, majorization characterizes when 
one quantum bipartite pure states can be transformed to 
another deterministically via local operations and classi- 
cal communication 20]. More details about majorization 
see |2y. 

In £| and 0, the authors related the majorization the- 
ory to quantum algorithms as follows: let \^ m ^) be the 
state of the register of a quantum computers at an op- 
erating stage labeled by m = 1, ...M, where M is the 
total number of steps in the algorithm. Let N be the 
dimension of the Hilbert space. Suppose {l?)}^! is the 
basis in which the final measurement is performed. Then 
suppose in this basis the state \i/j( m )) is 



|^ (m) ) = 



N 
i=l 



(11) 



If we measure |^( m )) in the basis {\i)}f =1 , the probability 
distribution associated to this state is 



P 



(m) = 



(rn) I (to) i2 



(m)\|2 



(12) 



where i = 1,2, ...,7V. If p( m ) -< p( m+1 ) for every to, we 
say this algorithm majorizes step by step. 

Especially, in Q and the authors applied majoriza- 
tion theory to study quantum adiabatic algorithms. By 
numerical simulation the authors observed that local 
quantum adiabatic search algorithm shows step-by-step 
majorization. For global adiabatic search algorithm, they 
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observed that though step-by-step majorization doesn't 
hold if the adiabatic inequality in Eq.(4) is tight, step- 
by-step majorization will appear as long as the evolution 
of the hamiltonian becomes slower and slower. Note that 
these two quantum adiabatic search algorithms both be- 
long to the big class of quantum adiabatic algorithms we 
have discussed above. 



III. STEP-BY-STEP MAJORIZATION OF 
GUIDE STATES 

In this work, we discuss majorization in this class 
of quantum adiabatic algorithms described by Eq.(7)- 
Eq.(lO). Firstly let's consider guide states of these al- 
gorithms. We prove that, for any quantum adiabatic 
algorithm of this big class, step-by-step majorization of 
the guide state holds perfectly. 

Theorem 1 Suppose H and Hi given by Eq. ( 7) and 

Eq. ( 8) are the initial and the final hamiltonians of a 
quantum adiabatic algorithm. Suppose this quantum adi- 
abatic algorithm has a linear path given by Eq.(9). Then 
the guide state of this algorithm shows perfect step-by- 
step majorization. 

Proof. Suppose the non-degenerate ground state of 

N 

H{s) = (1 - s)(I - \a)(a\) + *X)/(i)|i><i|, (13) 

i=l 



\4>(s)} = (ai,a 2 , ...,aAr) 2 



(14) 



and the corresponding eigenvalue is A(s). Suppose 
mini<j<Ar f(i) = 0. Otherwise we can let 

H(s) = H (s) - s X I X mm f(i), (15) 

which doesn't change the ground state of H(s). Without 
loss of generality, we suppose /(l) < /(2) < /(3) < 
... < f(N). Because if we change f(i) to f'(i) = f(n(i)), 
where it is a permutation of 1, 2, N, the ground state of 
H(s) for any s € [0, 1] will be changed to (&i, 62, &at) T , 
where bi = a^uy This can be proved as follows. Let 



H'(s) = {1 - s)H + sH[, 



where 



N 



(16) 



(17) 



i=l 



Note that P 7r if'(s)P 7r = H(s), where P w is a permutation 
matrix such that P^\i) — | tt (i)) . So we have 



Pr(&i, b 2 , b N ) T = (01, a 2 , a N ) J 



(18) 



That is to say, bi = . We will find that a permutation 
of the components of ground states doesn't affect our 
analysis later. 

By the definitions of \ip(s)) and A(s), we have 

ffOOIVOO) = A(*)M«)>. (i9) 

Substituting Eq.(13) and Eq.(14) into Eq.(19), we get 

N 



1 



^ ai = t(s)a l + sf(i)ai, 



(20) 



i=i 



where t(s) = 1 — s — A(s). Note that t(s) is the biggest 
eigenvalue of 



G(s) = (1 - s)I - H(s). 



(21) 



For every s S [0, 1) and ds > 0, an explicit calculation 
shows that 



N 



G(s + ds)-G(s) = -dsx(\a)(a\ + ) f(i)\i)(i\). (22) 



Because G(s) — G(s + ds) is a strictly positive matrix, 
t(s) is a strictly decreasing function of s 21] . It is easy 
to get t(0) = 1 and t(l) = 0. Then we have < t(s) < 1 
for any s £ (0, 1). 

By Eq.(20) we can obtain 

(t(s)+sf(i))ai = (t(s)+sf(j))a ] , i,j = 1,2,..., TV. (23) 

Thus it can be proved that if we let a\ > (Note that, 
cA(s) is also the ground state of H(s), where |c| = 1), 
we have > for any s € (0, 1) and any i. Otherwise 

if Oj = for any i we will have X)i=i a i 2 = accord- 
ing to Eq.(23), which is a contradiction. In this paper, 
we suppose the components of all ground states are real 
numbers. Substituting Eq.(23) into 



N 
i=l 



a, 



gives 



ai 



^AT / 1 y;' 



(24) 



(25) 



Note that t(s) is a strictly decreasing function of s, which 
means a\ is a strictly increasing function of s. However, 
for other ai, the monotony is a little more complicated. 
It's possible that they are not monotonous. However, 
we can prove that their increase and decrease are well- 
regulated. Namely, for s > and 1 — s > ds > 0, if 
o-i > a/, we have aj > a/, where (<xi', 02', a_/v') T is 
the ground state of H (s + ds) and i < j. 

Let's prove this conclusion. From Eq.(23), we get 



Oi = t(s) + sfjj) 
aj t(s) + sf(i) ' 



(26) 
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and 



Oi_ = t(s + ds) + (s + ds)f(j) 
a/ t(s + ds) + (s + ds)f(i)' 



(27) 



Because t(s) is a strictly decreasing function, it can be 
checked that 



*(*) + sfU) „ *(* + *0 + (« + <fe)/(j) 



< 



t(s) + s/(«) t(s + ds) + (s + ds)/(i) 



So 



a ; 



< 



(28) 



(29) 



Thus if di < cii, we have a,j' < aj. 

According to the discussion above, we know that for 
every s £ (0,1) there is a special integer io(s). When 
i 5~ *o( s ) we have dj < a/ and when i > io(s) we have 
di > a/. 

Now we are in a position to prove our first main con- 
clusion. Namely, 



(ar 2 " 2 



,a N 2 ) T -< {ai 2 ,a 2 ' 2 , 



, a N 



I2\T 



(30) 



Firstly, according to Eq.(23) it can be checked 



that the 
(ai , d2 , 



l«i 2 ,a2 2 , 



, a N 



components of 

.,ajv' 2 ) T are in decreasing order. 



2\T 



ondly, for any s £ 



(0, 1) and any k 

k k 



and 
Sec- 
1,2,...,N, if 



k < io(s): we have Yl a i — 2 a f because a, < a- 

i=l i=l 

for i < k. If fc > io(s), a, > for i > k, so we have 

AT JV k k 

E a 2 > £ af. Thus we also get E «f < E a? 



i = A:+l 



1. That completes the proof 



N 

because E a ? = E a 'i 

i=l i=X 

of Eq.(30), namely, step-by-step majorization holds if 
the adiabaticity is perfect. □ 
Note that if the form of Eq.(7) doesn't change, the 
ground state \a) in Eq.(7) can be replaced by any other 
vector of Hadamard basis. Because it can be proved if 
| a) is replaced by any other vector of Hadamard basis, 
for any s any component of the ground state of H(s) will 
not change up to the sign. Thus replacing the ground 
state in Eq.(7) doesn't destroy majorization. Besides, 
the path along which the hamiltonian varies in Eq.(9) 
can be replaced by any interpolation path in Eq.(2) pro- 



idcd 



is a increasing function of t, which doesn't 



/(t)+9(t) 

destroy step-by-step majorization either. 



IV. STEP-BY-STEP MAJORIZATION OF 
ACTUAL STATES 



state. We show that step-by-step majorization of the ac- 
tual state holds approximately 0, . 

Suppose in actual adiabatic evolutions, the actual state 
of the system is 



W(s)) = {b 1 M,-,b N ? 



Let 



B k = 



E 



k = 1,2, ...,7V. 



(31) 



(32) 



In |j] the authors studied s — B\ curve [B\ is the prob- 
ability of finding the right solution) and s — Bi curve of 
global quantum adiabatic evolution for search problem 
by numerical simulations. If step-by-step majorization 
holds perfectly, the curve should be monotonous. How- 
ever, they observed that oscillation appears at the end 
of s — B\ curve and s — Bi curve, which destroys step- 
by-step majorization (See Figure. 1). Furthermore, they 
also observed that the oscillation becomes weaker and 
weaker and step-by-step majorization appears as long as 
the running time becomes longer and longer. 




FIG. 1: A case that oscillation appears at the end of s — B\ 
curve (the solid curve). The dashed curve is s — A\ curve. 

Now, based on Theorem 1 we prove that for any 
quantum adiabatic evolution belonging to the big class 
discussed in this paper the oscillation at the end of 
s — Bk(l < k < N) curve, if any, will disappear if the 
running time becomes longer and longer. 

We consider an arbitrary state of the system near the 
end of the quantum adiabatic evolution. Let 



In the section above, we consider the guide state of a 
quantum adiabatic algorithm. In this section, we con- 
sider the actual state of the system. According to quan- 
tum adiabatic theorem, at any time of the whole com- 
putation process the actual state approaches the guide 



-4 t 



,. /■• 1.2 V. 



(33) 



where (oi, a./v) T is the guide state of the quantum 

adiabatic algorithm as before. Then from Eq.(23) and 
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Eq.(25) we have 

^ SrH 7 ^ , k = l,2,...,N, (34) 

2^=1 { l+f /(i) ) 

From Eq.(9) we get 



(35) 



1 - s " ' 1 - s 

It can be seen that yzj, the ground state eigenvalue of 
4j-ij, is a strictly decreasing function of s |2l|. So 



as 1 — s 



(36) 



which makes 



<i , t d , 1 — s — A . , 

= H < o. (37) 



g?s 1 — s ds 1 — s 

A simple calculation shows that 

dt t 
-— > 



ds 1 — s ' 



(38) 



then 



#(|)>t(- + T J -)' °< s<1 - ( 39 ) 
ds t ts 1 — s 

Calculating the derivative of Eq.(34) we get 

dA k * ^ /a|(f) /(o^(f) . 

i=2 j=fe+l 
iV 



+ f/(i) i + f/W 



,/(<)£(f) 



i=fe+l 



Since 



/0U(f) /«£(f) 



i + f m i + f /(*) 

when N > j > i > 1, we have 



ds 



N 



>2a\ £ <${ 



i=k+l 



+ im 



Let 



If m > 1, 



m = min f (i). 



/Wf(f) > f/«(7 + T i l) 



(41) 



(42) 



(43) 



1 + f/W f/(i) + l 
, gg + f /(*) 

i + f/W 



(44) 



>1. 



Here, we use Eq.(39) and the fact t < 1 < /(i) and 

ik >L . 

Similarly, if m < 1, 

/(0£(f) > f/W(i + T^) 



i + f m i /(*•) + 1 



i /(») i j_ 

=m- 



>m. 



So, we have 



/(0£(f) 
i + f/W 



> c, 



(46) 



where c = mzn{m, 1}. 

Substituting Eq.(46) into Eq.(42), we get 



dA k 
ds 



N 



>2c-a\ a i 

i=fc+l 



Note that 



We finally have 



a, > 



as re 



(47) 



(48) 



(49) 



According to quantum adiabatic theorem, we know 
that for any positive S we have a finite running time T 
such that 



\W(s)Ms))\>l-5 2 /2 



for any s £ (0, 1). 
Since 



m-w)\ 2 =2-2<vooi^(*)> 



it can be seen that for any s 

k 



(50) 



(51) 



E|a, ; -6 t | 2 <S 2 . 



(52) 



Here, we choose the global phase of \ip'(s)) such that 
(4>(s)\ip' (s)) is real. According to Cauchy's inequality, 



Because 



| a j — fej | < Vkd. 



\a,i + bi\ < 2, 



(53) 



(54) 
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we have 



k 

i=i 



b 2 i \=J2\a i -b i \-\a i + b i \<2VkS. (55) 



i=l 



Thus 



E n 2 - E n 2 i ^ E i« 2 - 52 i < 2 ^ 6 - ( 56 ) 



Now let us consider two points (si,^) an d (s2,A k ) 
on s — A k curve (about the guide state) and two points 
(si,B' k ) and (s2,B k ) on s — B k curve (about the actual 
state), where — S\ = As, < As <C 1. These four 
points are all near the end of the quantum adiabatic evo- 
lution. If step-by-step of the actual state holds, s — Bk 
curve should be a monotonically increasing curve. Sup- 
pose that Eq.(50) holds. According to Eq.(56) we have 



\A k -B k \ < 2Vkd, 



\A' k -B' k \ 



< 



2VkS. 



(57) 



On the other hand, according to the discussion above, we 
have 



Then if 



A k -A' k > ^A k (l-A k ).As, 



(58) 



(59) 



we have B' k < B k . Note that for arbitrary small As we 
can find corresponding 5 or running time T such that 
Eq.(59) holds. Thus, it can be judged that when the 
running time becomes longer and longer the oscillation 
at the end of s — B k curve becomes weaker and weaker, 
and at the same time step-by-step majorization appears 
slowly. This explains the results of numerical simulations 
in |4(, which is a special case of our discussion above. 
This is consistent with our intuition. Because of quan- 
tum adiabatic theorem, we know that when the running 
time becomes longer and longer, the distance between the 
actual state of the system and the guide state becomes 
smaller and smaller. Since it has been shown that the 
guide state shows perfect step-by-step majorization, step- 
by-step majorization of the actual state appears slowly 
when the adiabaticity of the quantum adiabatic evolution 
becomes better and better is natural. 



V. DISCUSSIONS AND CONCLUSIONS 

The results of this paper can help us to understand 
the relation between the efficiency of quantum algorithms 
and step-by-step majorization Usually, except some 



simple cases, a decisive mathematical analysis of a quan- 
tum adiabatic algorithm is not possible, and frequently 
even the estimation of the running time is very difficult. 
However, the minimal running time of the class of quan- 
tum adiabatic algorithms discussed above can be esti- 
mated, which is exponential in n, where n is the prob- 
lem size ^3, In this connection, the efficiency 
of these algorithms is not very good. However, though 
this class of quantum adiabatic algorithms can't offer 
exponential computational speed-up, they may be used 
to design quantum algorithms more powerful than cor- 
responding classical ones. In |9j, the authors resolved 
quantum search problem via a quantum adiabatic algo- 
rithm of this class, which is optimal and gives rise to a 
quadratic speed-up. Similarly, In also by a quantum 
adiabatic algorithm of this class the author tried to solve 
the hidden subgroup problem. It can be shown that the 
latter is not optimal. So it is clear that the efficiency of 
some of this class of algorithms is different from that of 
others. 

On the other hand, we have shown that if the run- 
ning time of these algorithms becomes longer and longer 
they all show step-by-step majorization. Thus this re- 
sult can be looked as a strong evidence that step-by-ste 
majorization is not a sufficient condition for efficiency 

In conclusion, we have shown that for any algorithm of 
a big class of quantum adiabatic algorithms, step-by-step 
majorization of the guide state holds perfectly. We also 
have shown that step-by-step majorization of he actual 
state holds approximately. The longer the running time 
is, the better majorization holds. This shows that ma- 
jorization seems to appear universally in quantum adia- 
batic algorithms. If the initial and the final hamiltonians 
are diagonalizable in Hadamard basis and computational 
basis respectively, does every quantum adiabatic algo- 
rithm show step-by-step majorization provided the run- 
ning time is long enough? We have shown that majoriza- 
tion is not sufficient for the efficiency of quantum algo- 
rithms. A natural question is, is step-by-step majoriza- 
tion necessary for efficiency? In £| the authors showed 
that any known efficient quantum algorithm does obey 
step-by-step majorization. These questions need further 
research. 
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